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ABSTRACT 

' The goal of this paper is to provide a numerically fast and stable description for the microlensing 

, amplification of an extended source (either uniform or limb-darkened) that holds in any amplifica- 

■ tion regime. We show that our method of evaluating the amplification can be implemented into a 

' light-curve fitting routine using the Levenberg-Marquardt algorithm. We compare the accuracy and 

computation times to previous methods that either work in the high-amplification regime only, or 
require special treatments due to the singularity of elliptic integrals. 

In addition, we also provide the equations including finite lens effects in microlensing light curves. 
We apply our methods to the MACHO-1995-BLG-30 and the OGLE-2003-BLG-262 events and obtain 
results consistent to former studies. We derive an upper limit for the OGLE-2003-BLG-262 event lens 
size. 

We conclude that our method allows to simultaneously search for point-source and finite-source 
microlensing events in future large area microlensing surveys in a fast manner. 
, Subject headings: daiYva'Attev — gravitational lensing — galaxies: halos — galaxies: individual (MS 1, 

CIh' NGC 224) — Galaxy: halo — galaxies: luminosity function, mass function 

o 

' 1. INTRODUCTION 

In large area microlensing surveys, one has to search for microlensing signatures in billions of variable sources. This is 
straightfor ward to do and co mputationally inexpensive in the point-source approximation. One either fits a Pa czyhski 
fSJ . light curve (|Paczvhskilll986f ). or, if appropriate, the Gould high-amplification approximation for point sources (jCouldl 
> ■ [1991 . One major disadvantage of these point-source light curves is the infinite amplification for a lens exactly in front 
^sO \ o f the point so urce. 

iGouldl (|1994l ) extended Paczyhski's light curve to finite sources which also avoids infinite amplifications. His equation 
CO ' describes the amplification as the two-dimensional integration of the Paczyhski amplification over the circular source, 
assumed to have constant surface brightness. 
. Using the limiting form of Paczynski light curve under high amplification, Gould is able to factor out the two - 
' dimensional integral into point source amplification times a much simpler integral. Meanwhile, IWitt fc Mad (jl994f ) 
^\ \ obtained the finite source amplification directly from the lens equation by comparing the area of the source and its 
i lensed images. However, one needs to take care of the singular points for the elliptic integrals of the first and the third 
kind when using their formula. 

In this paper we adopt the same strategy as iGouldl ()1994[ ) , because in this way more general surface-brightness 
profiles for the sources (e.g. limb-darkened ones) can be taken into account straightforwardly. For the cases of a 
uniform disk we will also compare our results with lWitt fc Maol ()1994l ). 
, The technical issue of the integration in the Gould extended source formalism can be carried out in several different 
ways. The two straightforward ones are to use polar c oordinates and to cho ose the coordinate ce nter either (1) at the 
source center or (2) at the lens center. IGouldl (|1994( l took the first choice. iRiffeser et al\ (j2006t ) have shown for the 
very special case where the lens is positioned along the line of sight to the source that the integration can be solved 
very easily if the second option is chosen. 

This leads us to choose the lens center as the coordinate center in general to benefit from the more simple integrand. 
We will show (in Section 2) that in this way the amplification of a uniform circular source is reduced to a one-dimensional 
integral and can be computed numerically fast and stable by using the composite Simpson's rule. A limb-darkened 
source is treated in Section 3. The two-dimensional integral can be solved numerically again in a fast and stable 
fashion, and light -curve fits for limb-darekened profiles can be obtained with the Levenberg-Marquardt algorithm (see 
iPress et aZ.ll2007f l with less than 100 ste ps. We also allow for finite lens sizes in Section 4. As a test example, we 
apply our fit ting methods to a M ACHO (jAlcock et al.lll99^ event and an OGLE (The Optical Gravitational Lensing 
Experiment: lUdalski et anil992f l event in Section 5. We conclude in Section 6. 

2. THE FINITE-SOURCE MICROLENSING EQUATION 

We first introduce our notation. Let Re be the Einstein radius of a point mass lens, and b be the impact parameter 
of a point source. Then one can write the amplification of the point source by the point mass lens as a function of the 
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Fig. 1. — Geometric definitions. Left: source is overlapping the lens center. Right: lens is outside the source radius. 

dimcnsionless impact parameter u := b/ Re = Q/Oe as 

^Ps(w) 



(1) 



u\/v?' + 4 
(jPaczvnskilll"986l ). 

If the source is extended one can obtain the lensed flux and the total amplification by integrating A(u) over the 
source area, weighted by the surface-brightness profile of the source. We now derive the amplification for a circular 
source^ with radius p {p= j^* is the projected source size in units of the Einstein radius Re and i?* is the physical 

source size). The situation is sketched in Figured] There arc two cases: either the center of the lens L (projected 
along the line of sight) is within the extended source centered at S (the left side of Figure [T|) or the lens is outside the 
extended source (the right side of Figure [T]), i.e., either u < p or u > p. One obtains the amplification of the extended 
uniform source by integrating the point-source amplification over the source area .Asourcc^ 
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Using polar coordinates centered on the lens, one can write 

IT "2(1?) 



A*{u-p) = ^^ J A^,{u)ududi} 

ui{^) 



The integration boundaries ui and U2 are 
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and so the amplification becomes 
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which can be approximated numerically using the composite Simpson's rule with 7T,(an even number) grids: 
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where = [112(1?) 1/1*2(1?)^ -I- 4 — ui(i?)\/iti(t?)^ + 4]. The upper limit of i? changes from -I to tt when the lens 
crosses the edge of the source from outside to inside, thus we set a grid of 2n for u < p in order to have the same step 
size on both sides. 

Gould (1994) argued that the finite-source effects are prominent only when the lens is very close to the source center 
(u ^ 1), and thus one can approximate Equation ([1]) by 



The reader is referred to lHevrovskv fc Loebl l|19971 ) for a more general case of elliptical source. 
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,ti < 1 



(8) 



and the finite-source light curve can be obtained by solving elliptic integrals (see also lYoo et aZ.|[2003 : ICassan et al\ 
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np 

where E{(j), k) is the elliptic integral of the second kind and -d 

i9 



is defined as 



f ^ ,u<p 
arcsin(/9/u) ,u>p 



(9) 



(10) 



We no w compare our method for A* {u; p) with previous ones, i.e. with lPaczvhiy (|T9M ). [GouO(|T99l and lWitt fc Mad 
()1994l ). and illustrate these comparisons in Figures [2] and [3] 

Equation ^ and (jlOp allow a fast computation of finite-source light curves in the Gould approximation, which however 
is accurate only for a high-amplification event. This is shown in F igure [H wher e for h igh amplifications (right panel) 
the Gould finite-source approximation {gray) is very close to the IWitt fc Maol ()1994[ ) light curve (displayed in solid 
black) ^ but fairly off when the lens transits the source for low amplifications (left panel). 



Moderate amplification regime, uo=0.1 



Higli amplification regime, uo=0.01 




-20 -15 -10 -5 



time [day] 

Fig. 2. — Comparison of finite-source light-curve approximations. Left: 
p = 0.5. Right: high-ampUfication regime with Ie = 10, uq = 0.01 and p 



time [day] 



moderate-amplification regime with Ie = 10, uq =0.1 and 
: 0.05. In dash ed bla ck the Paczyhski light curve for a point 
source, in solid black Witt & Mao light curve, in gray the approximation derived W' Gould ( 19S3J) and in dashed white Equation ((TJl with 
n = 10. The vertical lines indicate the time when u = p. Our formula is as good as [Gould ( 1994) in high-amplification regime and is better 
in the moderate-amplification regime. 



Our formalism from Equation ^ and that of I Witt fc Mad ()1994[ ) both provide the exact light curves for uniform 
extended sources. In the Witt fc Mao formalism one has to evaluate an elliptic integral which shows singularity when 
the impact parameter u is similar to the source size p. Witt fc Mao therefore derived a separate solution for the case 
oi u = p. This method is difficult to implement into numerical fitting routines in general, and particular cumbersome 
for those fast numerical fitting routines, where the partial derivatives have to be provided. 

We therefore suggest to start from our exact formalism given in Equation ^ and estimate values for the integral 
usi ng Equation Q wit h n = 10. The comparison with results from higher values for n = 500 or the comparison with 
the lWitt fc Mad ( 19941 ) formalism - see the gray and dash- dotted curves in Figure [3]- shows that Equation ([7]) (with 
n = 10) provides a precise numerical estimate for the integral already. Another advantage of our formalism is that 
one can obtain the derivatives of Equation ([5]) with respect to source radius p and m in a straightforw ard manner 
(see A ppendix A). This enables us to use fitting routines as, e.g., the Lcvenbcrg-Marquardt algorithm (see lPress et al\ 
[20071 ) which converge in t his cas e in les s than 100 iterations. 

The approximation by iGouldl (jl994f ) with A^^ (u) evaluated according to Equation ^ is actually valid for all u 
provided that p ^ 1, so it deviates from Equation ^ for larger source size. In fact, more than 80% (2548 out of 
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Fig. 3. — Percentage deviation in amplification compared to Witt &z Mao formalism The expression of IGouldl l|1994l') is valid for 

small source {solid black) but shows deviation > 2.5% for larger source {dotted black). Equation JTJl with n = 10 shows a smaller deviation 
(< 0.5%). Equation (0 with n = 500 for both source sizes are well overlapped with each other, so we show here only p = 0.05. 



3153) of the microlensing events detected from the OGLE experiment^ (jUdalskill2003f l have maximmn amphfication 
< 10 (see Figure 13]). This highUghts the necessity of a fast fitting routine for the moderate-amplification regime. We 
then compare the light-curve computation time of Equation (O to Gould's formalism (see Figure [5]). With n = 10, 
Equation ([7]) is about 38% faster then Gould's formalism when u < p and is > 55% faster when u > p. Therefore, 
our ri = 10 approximation turns out to be a practical fast fitting routine for both moderate- and high- amplification 
regimes. 



Amplification of OGLE events from 2002-2007 
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Fig. 4. — Maximum amplification of microlensing events detected 
by the OGLE experiment from 2002 to 2007. Most of the events 
(> 80%) have maximum amplification < 10. Events with maximum 
amplification > 100 , which are categorized into interval 100-110 in 
this plot, are relatively rare (< 4.4%). 
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Fig. 5. — Light-curve computation efficiency. We compare light- 
curve computation time of Equation jTj with n = 10 to that of 
Gould's formalism for various source radii (p = 0.01, 0.1, and 1). 
The computation time for our approximation is comparable to the 
Gould formalism; it is about 38% faster when u < p and is > 55% 
faster when u > p. 



3. FINITE SOURCE WITH LIMB DARKENING 

The next step towards a more precise microlensing light curve for extended sources is to account for limb darkening. 
Since the darkening is increasing towards the edges of the source, the limb darkening brings finite-source light curves 
closer to the Paczyhski light curve which can be considered as the most extreme limb-darkening model with a delta 
function. 

We use the one-parameter linear limb-darkening profile from lYoo et al\ ()2004l ) for the surface brightness of the 
source, 



5(r/i?*,r) = S 




(11) 



^ http:/ /ogle. astrouw.edu.pl/ogle3/ews/cws. html 
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where r is the distance to the source center. 

r is the limb-darkening coefficient, and depends on the wavelength range used for the observations. S is the mean 
surface brightness of the source and defined as 

27r R, 

S{r/R^,T)rdfdd ^ nR^'^S . (12) 





We implemented the limb-darkening effects in our finite-source light curve as follows: 
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(13) 
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Equation p3p is still a double integral over u and 1?. But even here the divergent part cancels, and the function is 
numerically stable and can be evaluated using a small grid. The limb-darkening effects under moderate-amplification 
regime is shown in Figure [51 
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Fig. 6. — Limb-darkening effects on the finite-source light curve in the modcratc-amphfication regime. In dotted black we show the 
Paczyhski light curve for a point source with t^ = fO and tig = 0.1. In solid black, we show the finite-source light curve for a uniform 
source with a projected source size of p = 0.5. In dashed line and dash-dotted line, we plot the limb-darkened finite-source light curves with 
r = 0.3 and 0.6. Increasing V enhances the limb-darkening effects thus brings the finite-source light curve closer to Paczyrisky's formalism. 



4. FINITE-SOURCE EQUATION WITH FINITE LENS 

Given a finite-size lens, one can always find a time interval when the lens obscures the inner (and the outer, depending 
on the lens size) lensed image in the early rising stage and in the final declining stage of the light curve. In the following, 
we investigate how large this effect is depending on the lens size. 

Agol (200^ derived the lens-modified amplification by calculating how much area is unobscured by the lens in the 
image plane. One has to solve for the image position by inverting the lens equation and one has to evaluate the image 
area from the image boundary using Stokes' theorem. Depending on the source and lens radii, there are 7 different 
cases for the inner image and 6 cases for the outer image to be considered if one follows the derivation of Agol. 

Here we show that the finite lens amplification of a finite source again can be much more easily evaluated if one uses 
the polar coordinates u and d again. First, we consider a lens with physical radius i?ions transiting the surface of the 
source. The light emitted at a given point from the source follows the lens equation 

_^^h^_RE_ 

Re Re b± 

which gives the position of the two images in the lens plane (recall u 

b± _ u± \/v? + 4 
'R^ " 2 
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Re 



(14) 
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b_ obscured by lens 



Fig. 7. — Image obscuration by a finite lens with radius Riens = 0.5Re- 



b. obscured by lens 



with amplifications 
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2M^^2 
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= ^2 



(16) 



Here, denotes the outer image, and denotes the inner image in units of the Einstein radius. The sum of 



A+(u) and A-(u) gives the Paczyhski hght curve. 



An image is unobscured if —4— > Pions or > pions holds, where p\c 



pions fl.^ ^ Anions ...^v..^ pions — j^^ 

Einstein radius. Following this critcrium and Figure[7l there exists an upper limit for < 



is the lens radius in units of the 



B.E 



-pic 



and a lower limit for 



> Picns to be unobscured by the lens. Therefore, we only need to consider these two limitations when integrating 



the amplification in Equation ([3]): 



u du d-d 







[U1(<?),P1C.„3-^] 2 



min[M2(i?),-Picns + 5^] 



(17) 



d§ 



where Q{x) defines the Heaviside step function. 

Combining Equation (fT3|) and Equation (fT7|) fully considers a limb-darkened source and a finite lens: 



p) = ^ j j 



A+{u)Q ( - pions") + A_(M)e f-;^ - pic 



S'(r/i?,,r)udud79 . (18) 







5. RESULTS 

When we implemented the finite-source fitting using the Lcvenberg-Marquardt algorithm, we recognized that a good 
set of initial values is needed to bring the algorithm to convergence. Fitting a Paczyiiski light curve to derive these 
initial values for the finite-source fitting leads to very good results. The algorithm is stable and for an initial value of 
p = 0.1 it conv e rges w ithin 100 iterations. 

lAlcock et ai\ (fl99l were able to measure a microlensing light curve with finite-source effects in MACHO-1995-BLG- 
30. We extracted the data points from their paper and applied our finite-source fitting algorithms to them. Fitting 
Equation (|7l) with rt = 10 to the data yields a perfect agreement (see Tabic [T] and Figure [8]) with the parameters given 
i n Table 2 of lAlcock al\ (|1997f ): 

lAlcock et al\ ( 19971 ) t hen obtained the limb darkening coefficients of MACHO-1995-BLG-30 utilizing spectroscopic 
information. However, iHevrovsk?! ()2003f ) argued that the surface-brightness profile of this event can not be fully 
recovered due to its intrinsic complex variability. Therefore, we tested our limb-darkening fitting routine to another 
limb-darkened fi nite-source event OGLE-2003-BLG-262. Our results are shown in Tabled Figure [H] and Figure [TU] in 
comparison with lYoo et'ai\ (|2004[ ). 
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Fig. 8. — Finite-source light-curve fits for MACHO-1995-BLG-30 assuming a uniform source. Data points in R are from MACHO {red), 
CTIO, UTSO, WISE, and MJUO (gray) and V are from MACHO (blue) and UTSO (green). The dashed line shows the fight curve for a 
point-source model. The best-fitting finite-source light-curve parameters are displayed in Table [T] 



TABLE 1 

Light-curve parameters for MACHO- 1995-BLG-30 



Fit 


Ap^ (u) (this work) 


A*{u;p) (this work) 


Alcock et al. (1997), Table 2=* 


to 


1321.260 ± 0.002 


1321.235 ± 0.002 


1321.2(1) 




34.41 ± 0.02 


34.25 ± 0.02 


33.68(1) 


uo 


0.04133 ± 0.00004 


0.05569 ± 0.00006 


0.05579(1) 


P 




0.0722 ± 0.0001 


0.07335(1) 



OGLE-2003-BLG-262 fitting residual 
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Fig. 9. — Residuals of the observed light curve relative to the best-fitted point-source light curve. The solid black curve shows the light 
curve of an extended source with uniform surface brightness. The solid blue, solid red, and solid green curves are extended source models 
incorporating limb darkening in V, I and H bands with (Ty ,r i ,r u) = (0.72,0.44,0.26). The vertical lines indicate to for the best-fitted 
point-source (solid) and limb-darkened finite-source (dashed) model. For the light curves with the limb-darkened source we have left to as 
a free parameter. The best-fitting value for to slightly differs (see Table [2)l. This causes the asymmetric pattern of the residual relative to 
the Paczyhski light curve. 

Finally, we choose several lens sizes for the configm-ation of OGLE-2003-BLG-262 to investigate the influence of the 
finite lens effects on the microlcnsing light curve in Figure [101 The light curve is strongly altered only if the lens size 
is comparable to or larger than the Einstein radius. 

When the lens size is smaller than 0.93i?£;, it only partially covers the outer image and the finite lens effects can be 
observed only at the very beginning of the rising and near the end of the declining stage of the lensing event. Therefore, 
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TABLE 2 

Light-curve parameters for OGLE-2003-BLG-262. 



Fit (u) 



to 

p 



2839.852 ± 0.001 
12.83 ± 0.01 
0.02877 ± 0.00008 



A*{u;p) 

2839.838 ± 0.001 
12.61 ± 0.01 
0.0365 ± 0.0002 
0.0581 ± 0.0002 



^*r,n(";P) 

2839.8361 ± 0.001 
12.559 ± 0.016 
0.0361 ± 0.0002 
0.0598 ± 0.0002 



Note. We fixed the limb-darkening coefficients at {Tv,Ti,rH) = (0.72,0.44,0.26) 
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Fig. 10. — Finite-source and finite-lens light-curve fits for 
OGLE-2003-BLG-262. Data points are in /(square), l/(circle), and 
(triangle). The dashed line shows the light curve for a point source. 
The solid line shows the light curve for an extended source with uni- 
form surface brightness. The dotted lines illustrate the effects of finite 
lens sizes on top of finite-source size for lens sizes of piens = 0.93, 0.96, 
0.99, and 1.02. 



HJD' 

Fig. 11. — Fitting residual of various lens radius relative to the 
point-source model. Data points are in /(square), l/(circle), and 
//(triangle). The solid black line shows the light curve for an extended 
source with uniform surface brightness. The dotted lines illustrate the 
effects of finite lens sizes on top of finite source size for lens sizes of 
Picns = 0.93, 0.96, 0.99, and 1.02. One sees that all these cases can 
be safely excluded. 



we fitted various lens sizes up to picns = I-^Re using the full OGLE /-band data set. However, no improvement in 
has been found by introducing lens sizes as an extra parameter in the finite-source model (see Figure [T^ . This implies 
that the lens size effect is negligible for OGLE-2003-BLG-262. 

6. CONCLUSION 

We have demonstrated that finite-source effects can be more conveniently evaluated in the Icns-centcrcd polar 
coordinate system. The uniform source case can be reduced to a one-dimensional integral, which can be solved in 
a fast and numerically stable manner. The previously available formalisms were either comparably fast but held 
only in the high-amplification regime (the Gould finite-source approximation) or held in any amplification regime but 
involved an integral which has singularity and is slower to solve (the Witt & Mao approach). We also showed that 
the vast majority of the OGLE-lensing events have maximum amplifications smaller than 10, and therefore cannot be 
precisely described in the high-amplification, finite-source approximation of Gould. Our formalism allows a fast and 
simultaneous search for microlensing events with extended or pointlike sources in any amplification regime. 

We also presented the limb-darkening effects and finite lens size effects in our formalism. We showed for the case of 
OGLE-2003-BLG-262 how one can constrain the source size and obtain upper limits for the lens size. 

The Appendix provides the partial derivatives of the amplification for a uniform surface brightness source (Appendix 
A), a limb-darkened source (Appendix B), and a uniform surface brightness source with a finite lens (Appendix C), 
which are required in, e.g., the Levenberg-Marquardt algorithm to obtain microlensing light-curve fits. 



We thank Johannes Koppenhofer for fruitful discussions. This work was supported by the DFG cluster of excellence 
'Origin and Structure of the Universe' (www.universe-cluster.de). 
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Fig. 12.— contour map of OGLE-2003-BLG-262. In white contour, levels for 1, 2, and Str for source and lens size fitting are shown. 
The black triangle indicates the best-fitted model with p = 0.06056 and -Riens = 0- This suggests that the point lens assumption is sufRcicnt 
for OGLE-2003-BLG-262. 
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APPENDIX 

A. PARTIAL DERIVATIVES OF THE FINITE-SOURCE AMPLIFICATION FOR A SOURCE WITH UNIFORM SURFACE 

BRIGHTNESS 
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f2 



u sin 'd 
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cosd 



u sin 1? 



(Al) 



B. PARTIAL DERIVATIVES OF THE FINITE-SOURCE AMPLIFICATION FOR A SOURCE WITH LIMB DARKENING 
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TT "2 , 

2 f f 11+2 



2 f f ^ 



TP" J J \/P+4 dTx 

Ml 



(It) 



cZm dd 



with^ = cos sin^ d/^p^ - sin^ i?, ^ = cosiJ-Msin^ ^/^p2 _ ^2 31^2^^ = -p/^p2 _ ^2 sij^2 



P/ V — sin^ 



^' (It) = -irA(-27icost9 + 2u)/ (f^T^ 



-2uu COS "d-yu'^ 



■u?'—2uu COS 

p2 



)' i(ft) = |rA(^22-2^.Scos^ 
when wi and U2 are not equal to zero. 



C. PARTIAL DERIVATIVES OF THE FINITE-SOURCE AND FINITE-LENS AMPLIFICATION ASSUMING A SOURCE WITH 

UNIFORM BRIGHTNESS 



OA" 



— (u,p,pions) 



^(U,p,pie„s) 



^(U,p,pie„s) 



7^/ 4^[(^ps("2) + i)6'("^-P><=- + 7^) + (^ps("^)-1)Q(-"2-Pio„= + 7;^)]"2 

-V;^ I ^ [('4ps("l) + l)e("l-Plon„ + ^)+(Ap3(Mi)-l)e(-Ml-p,_ + ^)]Mi d€ 
7^/ ^[('4ps("2) + l)e(M2-Plc„. + ^)+(Ap3(M3)-l)e(-M,-pi_-|-^)]M, 

I ^ [(^Ps("i)+i)e("i-P><=-+7;ii7)+(^Ps("i)-i)Q(-"i-pio„=+^)]"i d.5 



(CI) 



1 -^ (".P.Plcns) 



+ ^/ /(^Ps(")-l)<-"-P.™=+^) 



- Plena + 



8Plo 



^'""^ ^ M dM diJ 



The derivatives can be obtained numerically by utilizing the same approache as shown in Equation ([7]). We also find 
that for u > p, substituting integration variable with v = ^ sin gives a numerically more stable estimations of the 
derivatives for a uniform brightness source: 

1 



dA* 





1 



(M''-p^)(j7+a)-4a _ (M^-p^)(c/-n)+4a 



dv 



with U = \/ — w2p2 and il = 1 — 



^((7+0)2+4 ^(!7-J2)2+4 

M^(M^-p^)(c/+r!)-4(OM^p^ M^(u^-p^)(E/-o)-i-4(nM^p^+;7-') 



^(!7+0)2+4 



V(C/-n)^+4 



(C2) 
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